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Abstract 

We analyze the spin-orbit coupling term in the nuclear energy density functional in 
terms of a zero-range NN-contact interaction and finite-range contributions from two-pion 
exchange. We show that the strength of the spin-orbit contact interaction as extracted 
from high-precision nucleon-nucleon potentials is in perfect agreement with that of phe- 
nomenological Skyrme forces employed in non-relativistic nuclear structure calculations. 
Additional long-range contributions from chiral two-pion exchange turn out to be relatively 
small. These explicitly density-dependent contributions reduce the ratio of the isovector 
to the isoscalar spin-orbit strength significantly below the Skyrme value 1/3. We perform 
a similar analysis for the strength function of the (V/o) 2 -term and find values not far from 
those of phenomenological Skyrme parameterizations. 
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The microscopic understanding the dynamical origin of the strong nuclear spin-orbit force 
O ■ is still one of the key problems in nuclear physics. The analogy with the spin-orbit interaction 
in atomic physics gave the hint that it could be a relativistic effect. This idea has lead to the 
construction of the scalar- vector mean-field models for nuclear structure calculations [1, 2]. In 
these models the nucleus is described as a collection of independent Dirac quasi-particles moving 
in self-consistently generated scalar and vector mean-fields. The footprints of relativity become 
visible through the large nuclear spin-orbit coupling which emerges in that framework naturally 
from the interplay of two strong and counteracting (scalar and vector) mean-fields. The corre- 
sponding many-body calculations are usually carried out in the Hartree approximation, ignoring 
exchange terms as well as the negative-energy Dirac-sea. The nucleon-nucleon interaction intro- 
duced into these models is to be considered as an effective one that is tailored to properties of 
finite nuclei but not constrained by the observables of free NN-scattering. For recent ideas about 
how the large Lorentz scalar and vector mean-fields can be linked to the condensate structure 
of the QCD- vacuum, see ref. [3]. 

On the other hand phenomenological Skyrme forces have been (and still are) used extensively 
for the non-relativistic description of nuclei. In that approach the nuclear spin-orbit coupling is 
introduced by hand through a suitable NN-contact interaction (i.e. a spin-orbit delta-force) with 
an adjustable strength parameter Wq. Although there is a wide spread in the (other) parameters 
of the about hundred available Skyrme forces the one related to the spin-orbit coupling turns 
out to be rather stable with a value around Wq ~ 120 MeVfm 5 . In fact one may even ignore the 
concept of an underlying effective (zero-range) two-body interaction and just take the emerging 
(parameterized) energy density functional as the true starting point for nuclear structure cal- 
culations within the self-consistent mean-field approximation. Such an interpretation is then in 
line with the energy functional approach to many-body problems of Kohn and Sham. 

The more basic Dirac-Brueckner approach of refs.[4, 5] solves a relativistically improved 
Bethe-Goldstone equation with a realistic nucleon-nucleon interaction parameterized in terms 
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of one-boson exchange potentials. This Lorentz-covariant approach is able to describe proper- 
ties of nuclear matter and finite nuclei without any new adjustable parameter. In particular, 
the spin-orbit splittings of medium and heavy nuclei are correctly reproduced (in local density 
approximation). A drawback of the method is that the underlying NN-interaction has been con- 
structed without the explicit use of chiral symmetry. In a relativistic Hartree-description there 
is a close connection between the mean-field potential and the spin-orbit one, since they differ 
only by the sign of the vector-meson exchange contribution, while in non-relativistic approaches 
there is no connection at all. Unfortunately there is no direct experimental hint for either possi- 
bility. Despite all its known successes the fully relativistic treatment of the nuclear many-body 
problem seems unnatural in view of the small ratio (kf /M) 2 ~ 0.08, with kf ~ 263 MeV the 
Fermi momentum of equilibrated nuclear matter and M = 939 MeV the (free) nucleon mass. 

The purpose of the present paper is to show that the spin-orbit interaction relevant for the 
structure of (medium and heavy) nuclei is quantitatively consistent with that relevant for low- 
energy elastic NN-scattering. We demonstrate that the short-range spin-orbit interaction as 
extracted from various high-precision nucleon-nucleon potentials (which accurately fit all NN- 
phase shifts and mixing angles below the NN7r-threshold) is in perfect agreement with the one 
needed in non-relativistic nuclear structure calculations (using the Skyrme phenomenology). 
Additional long-range spin-orbit couplings generated by chiral two-pion exchange turn out to 
be relatively small. These explicitly density dependent effects reduce favorably the ratio of the 
isovector to the isoscalar spin-orbit coupling strength. We perform also a similar analysis of the 
(Vp) 2 -term and find reasonable values for its strength function (at densities p > 0.05 fm -3 ). 

Let us begin with writing down the explicit form of the spin-orbit coupling term in the nuclear 
energy density functional: 

£so[Pp, Pn, J P , Jn] = F so (k f ) Vp • J + G so (k f ) Vp„ • J v , (1) 

where the sums p = p p + p n , J — J p + J n and differences p v — p p — p n , J v — J p — J n of proton 
and neutron quantities have been introduced. 

P P ,n(r) = = E , (2) 

" aGocc 

denote the local proton and neutron densities which we have rewritten in terms of the corre- 
sponding (local) proton and neutron Fermi-momenta k p>n (f) and expressed as sums over the 
occupied single-particle orbitals ^ p a l(r). The spin-orbit densities of the protons and neutrons 
are defined similarly: 

J P ,n(r) = J2 ^ P f(r)i * x . (3) 

aGocc 

Furthermore, F so (kf) and G so (kf) in eq.(l) denote the density dependent isoscalar and isovector 
spin-orbit strength functions. In Skyrme parameterizations [6] these are density-independent 
constants, F so (kf) = 3G so (kf) = 3W /A ~ 90MeVfm 5 , determined by one single spin-orbit force 
parameter Wo- 

The starting point for the construction of an explicit nuclear energy density functional 
£so [Pp, Pn, J P , Jn] is the bilocal density- matrix as given by a sum over the occupied energy eigen- 
functions: J2aeocc ^ p a n (^~ ^/2)^ r p a ^(^'+ a/2). According to Negele and Vautherin [7] it can be 
expanded in relative and center-of-mass coordinates, a and r, with expansion coefficients deter- 
mined by purely local quantities (nucleon density, kinetic energy density and spin-orbit density). 
As outlined in section 2 of ref. [8] the Fourier-transform of the (so expanded) density-matrix 
defines in momentum-space a medium-insertion for the inhomogeneous many-nucleon system. 
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It is straightforward to generalize this construction to the isospin-asymmetric situation of dif- 
ferent proton and neutron local densities, p p , n (r ) and J Pjn {f). We display here only that part of 
the medium- insertion T(p, q ) which is actually relevant for the diagrammatic calculation of the 
isoscalar and isovector spin-orbit terms defined in eq.(l): 

r(p,q) = J dhe-^[^^e{k p -\p\) + l -^e{k n -\p\) 

'S(k f - \p\) - k f S'(k f - \p\)} (axp). (J + r 3 X)} . (4) 
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The double line in Fig. 1 symbolizes this medium insertion together with the assignment of the 
out- and in-going nucleon momenta p ± q/2. The momentum transfer q is provided by the 
Fourier components of the inhomogeneous matter distributions p Pt n{r) and J pn (r). 

Next, we write down the lowest order four- nucleon contact-coupling which generates a spin- 
orbit interaction: 

d$tr = -^{(JVtViV) • (VAft x ?iV) + (V^iV) ■ (iV 1 ? x W) 

—(N^N) (VN^a x ViV) + (N*BN) ■ (V7V f x ViV)} . (5) 

This Lagrangian, together with the notation C5 for the coupling constant, has been taken over 
from ref. [9]. The somewhat lengthy expression in eq.(5) is dictated by Galilei invariance as will 
become clear immediately. The contribution of the contact-Lagrangian to the T-matrix of 
the scattering process N(p 1 ) + N(p 2 ) — > N(p{) + N(p{) reads: 

C C 
-iy (^1 + ^2) • [{pi ~ P2) x (pi - p 2 )\ = (^i + <? 2 ) • [{pi - Pi) x (p\ - p 2 )\ , (6) 

where Galilei invariance is now manifest, since only differences of momenta occur. 0^2 denote the 
conventional spin-operators of the two nucleons. Furthermore, by comparing with the analogous 
T-matrix element of the Skyrme spin-orbit interaction proportional to Wo (see eq.(4.105) in 
ref. [10]) one finds the relation C 5 = 21^0 between coupling parameters. As a side remark we 
note that the spin-orbit term which arises from heavy scalar and vector boson exchange between 
nucleons has the property of Galilei invariance only if the respective ratios of mass to coupling 
constant are equal: ms/gs = mv/gv- The corresponding contribution to the strength parameter 
C 5 reads in that case: c!f +V ^ = 2(1 + K, v ){gv/Mm v ) 2 where M = 939 MeV stands for (average) 
the nucleon mass and Ky denotes the tensor-to-vector coupling ratio of the heavy vector boson 
(e.g. the a;(782)-meson). 

In the work of Epelbaum et al. [9] numerical values are given for the so-called spectroscopic 
low-energy constants which characterize the short-range part of the nucleon-nucleon potential 
in certain low partial waves (S- and P-waves). The spin-orbit strength parameter C 5 = 2W 
of interest here is determined by the following linear combination of the 3 P-wave low-energy 
constants: 

° 5 = itb [ 2C(3Po) + 3C(3Pi) ~ 5C ^} ■ (7) 

This simple relation is central to our work since it allows to connect the phenomenological 
Skyrme parameter W = C 5 /2 with the low-energy dynamics of elastic NN-scattering. 
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Fig. 1: Left: The double line symbolizes the medium insertion T(p,q) defined by eq.(4)- Next 
shown are generic two-body Hartree and Fock diagrams. Their combinatoric factor is 1/2. The 
wiggly line symbolizes the spin-orbit NN-interaction. 

We continue with the general form of the spin-orbit part in the NN T-matrix (using the 
sign-convention of ref.fll]). For the scattering process N(p) + N(—p) — > N(p + q) + N(—p — q) 
in the center-of-mass frame it reads: 



NN 



V so (q) + fi • r 2 W so {q) i{S x + <7 2 ) -(qxp), 



(8) 



Here, q is the momentum transfer between both nucleons and V so (q) and W so (q) denote the 
isoscalar and isovector spin-orbit NN-amplitudes, respectively. With the help of the medium 
insertion T(p,q) in eq.(4) one can now calculate diagrammatically the nuclear energy density 
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functional £ so [p p , p n , J p , J n }. From the Hartree and Fock diagrams in Fig. 1 one obtains the 
following formulas for two-body contributions to the isoscalar and isovector spin-orbit strength 
functions: 



FUkf) 



(2-body) _ 



= -i{3\/ so (0) + V so (2k f ) + 3W so (2k f ) + £dxx 



V so (2xk f ) +3W so (2xk f )]\ , (9) 
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H/ so (2A; / )-V; o (2A; / )-3W so (0) 



W so {2xk f ) - V so (2xk } 



(10) 
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Here, the terms — V r so (0)/2 and — W so (0)/2 belong to the Hartree diagram (with two closed nu- 
cleon lines) while the remaining ones summarize the contribution from the Fock diagram (having 
just one closed nucleon line). After the obvious identification V so (qY ct ^ = — C 5 /2 (compare eq.(6) 
with eq.(8)) one derives the following contribution of the four- nucleon contact- vertex to the spin- 
orbit strength functions: 



F so (k f )W = 
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G so (k f )W = 



The fixed ratio of isovector-to-isoscalar spin-orbit strength G so (kf)/F so (kf) = 1/3 is a conse- 
quence of the Pauli exclusion principle for a zero-range spin-orbit interaction. 

The long-range contributions to the spin-orbit NN-potential arise naturally from two-pion ex- 
change between nucleons in the form of a relativistic 1/M-correction. The corresponding lowest 
order 27r-exchange triangle and box diagrams have been evaluated in section 4.2 of ref.[ll]. In- 
serting the corresponding analytical expressions for V so (q) and W so (q) (see eqs. (22,23) in ref.[ll]) 
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into the master formulas eqs.(9,10) one gets the following contributions from chiral 27r-exchange 
to the spin-orbit strength functions: 



g 2 A m * 
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arctan u 
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G so (k f ) 



(27T) 



g 2 A m * 
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2m 2 



ln(l + u 2 ) 



+ 



u 



(2-5g A )+8u(l-4g A ) 



arctan u 



(13) 



with the abbreviation u = kf/m n where m n = 135 Me V stands for the (neutral) pion mass. As 
usual fa = 92.4 MeV denotes the weak pion decay constant and we choose the value g A = 1.3 for 
the nucleon axial vector coupling constant. Note that we have normalized the density dependent 
expressions in eqs.(12,13) to the value zero at zero density {kf = 0). This way one eliminates au- 
tomatically all (regularization dependent) short-range contributions which conceptually belong 
to the low-energy constant C 5 . As stressed in ref.[9] such a separation of long and short-distance 
dynamics is necessary in order to be able to compare with results extracted from realistic NN- 
potentials. At the next order in the small momentum expansion there are 27r-exchange diagrams 
with one chiral 7r7riViV-contact vertex. The corresponding spin-orbit NN-amplitudes have been 
written down in eqs.(13,15,16) of ref.[12]. After performing the necessary integration and sub- 
traction at kf — one gets the following 27r-exchange contributions to the spin-orbit strength 
functions: 



Fso(k f ) 



mt 



20u z 



c 4 + g 2 A (5c 4 - 2c 2 )J ln-^ 



3tt 2 M(4/ w ) 4 

— 29c4 + g\(58c2 — 6IC4) + u 2 gi(2c 2 — 5c 4 ) — c 4 



+ 



2 

-(l3c 4 + g 2 A (29c 4 - 26c 2 )) + 20w(c 4 + g 2 A (5c 4 - 2c 2 )) 



V1 + 



1// 



x ln(u + Vl + u 2 ) + — c 4 + g\[c 4 - 2c 2 ) ln 2 (u + Vl + u 2 ) , (14) 



Gso(kf) 



mz 



20w z 



C4 + g A ( 5c 4 + 6c 2y 



111 X 



9tt 2 M(4^) 4 

+29c 4 + ^l(174c 2 + 61c 4 ) + u 2 [c 4 + g 2 A (6c 2 + 5c 4 ) 
- -(l3c 4 + g\(78c 2 + 29c 4 )) + 20w(c 4 + g 2 A (5c 4 + 6c 2 )) 



VTT 



u- 



x \n(u + Vl + u 2 ) - — c 4 + g\{c 4 + 6c 2 ) ln 2 (u + Vl + u 2 ) . (15) 



The two occurring low-energy constants c 2 = 3.2 GeV -1 and c 4 = 3.4 GeV -1 [13] represent 
mainly effects from (single) virtual A(1232)-isobar excitation. In these expressions we have also 
introduced A; 2 -terms proportional to the chiral logarithm h^m^/A) in order to guarantee finite 
results for F so (kfY 2n ~ Cj ^ and G so (kf)^ 2n ~ Cj ^ in the strict chiral limit = 0. 

Now we can turn to numerical results. In Table I we have collected values for the isoscalar 
spin-orbit strength 31^0/4 of various phenomenological Skyrme forces, SIII [14], SkM [15], SkP 
[16], Sly [17], MSk [18] and Ski [19]. In the cases of Sly, MSk and Ski we have performed averages 
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Skyrme force 


SIII SkM SkP Sly4-7 MSkl-6 SkIl-5 


3W /4 [MeVfm 5 ] 


90.0 97.5 75.0 93.2 87.6 92.7 



Table I: Numerical values of the spin-orbit strength parameter 3Wo/4 for various phenomeno- 
l oaicaL Skyxme. faeces 



NN-potential 



Bonn-B CD-Bonn* Nijm-93 Nijm-I Nijm-IF AV-18* NNLO 



3C 5 /8 [MeVfm 5 



80.3 89.6 79.9 82.4 87.7 88.9 73 • • • 92 



Table II: Numerical values of the short-range spin-orbit strength parameter 3Cs/8 for various 
nucleon-nucleon potentials. The so-called high-precision potentials are marked by an asterisk. 

over several (slightly different) parameter sets Sly4-7, MSkl-6 and SkIl-5. One sees that the 
entries in Table I scatter somewhat around a mean value of about 3W /4 ~ 90 MeVfm 5 which 
can therefore be considered as the empirical spin-orbit strength relevant for nuclear structure. 
In Table II we give the equivalent spin-orbit strength parameter 3C 5 /8 as extracted from various 
realistic nucleon-nucleon potentials. These numbers have been produced with the help of Table 
IV in ref.[9] by forming the appropriate linear combination eq.(7) of C( 3 Po,i,2)- The last entry 
"NNLO" in Table II corresponds to the chiral NN-potential of Epelbaum et al. [9] where the 
cut-off A has been varied between 0.5 and 0.6 GeV (see Table I in ref.[9]). It is astonishing 
to observe that these realistic NN-potentials which are essentially adjusted to empirical low- 
energy NN-scattering data give numbers close to those of the phenomenological Skyrme forces. 
This holds in particular for the so-called high-precision NN-potentials CD-Bonn [20], Nijm-II 
[21] and AV-18 [22] (marked by an asterisk in Table II) where 3C 5 /8 comes out very close to 
3C 5 /8 ~ 90 MeVfm 5 . Therefore one can conclude that the strength of the spin-orbit interaction 
necessary for nuclear structure is perfectly consistent with that needed to describe low-energy 
NN-scattering data. This observation is the one of the main results of the present work. 

The recently constructed universal low-momentum nucleon-nucleon potential Vi ow _k offers 
another possibility to extract the short-range spin-orbit parameter 3C5/8. From the curvatures 
of the (diagonal) 3 P-wave potentials shown in Figs. 4,18 of ref. [23] we get after multiplying the 
relevant linear combination 2( 3 P ) + 3( 3 Pi) — 5( 3 P 2 ) with the conversion factor 37r/8M the value 
3C5/8 = 79 MeVfm 5 from Viow-k- This somewhat smaller value (in comparison to the high- 
precision potentials) finds its explanation in the fact that the "full" NN-potential and Vi ow -k 
differ by some local contact-terms [24]. Indeed the numerical values of the counter-terms listed 
in Table 1 of ref. [24] give rise to an additional a contribution to 3C 5 /8 of 10.4 MeVfm 5 which 
then closes the gap the empirical value of 3C 5 /8 ~ 90 MeVfm 5 . 
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Fig. 2: The isoscalar and isovector spin-orbit strength functions F so (kf) and G so (kf) versus the 
nucleon density p = 2kj:/3n 2 . The dashed curves show the result of the spin- orbit contact- coupling 
C5 and irreducible two-pion exchange eqs. (12,13). The full curves include in addition the two- 
pion exchange contributions eqs. (14,15) proportional to c 2i4 . The scale in the chiral logarithm is 
chosen as A = 0.5 GeV. 



In Fig. 2 we show the isoscalar and isovector spin-orbit strength functions F so (kf) and G so (kf) 
as a function of the nucleon density p = 2kj/37i 2 . The dashed lines result from the average value 
3C 5 /8 = 88.7MeVfm 5 of the three high-precision NN-potentials and the lowest order irreducible 
27r-exchange contributions eqs. (12,13). The full lines include in addition the next-to-leading 
order 27r-exchange contributions proportional to the low-energy constants c 2 = 3.2 GeV" 1 and 
C4 = 3.4 GeV" 1 [13]. The scale A in the chiral logarithm \n(m n /A) has been chosen equal to the 
momentum cut-off A = 0.5 GeV of ref. [9]. In the case of the isoscalar spin-orbit strength function 
Fso(kf) one finds that the two-pion exchange effects are relatively small and that they also cancel 
each other to a large extent. In contrast to this behavior does the isovector spin-orbit strength 
function G so (kf) get reduced with increasing density. At nuclear matter saturation density 
Po ~ 0.16fm~ 3 the ratio G so (kfo) / F so (kfo) is about 1/8, substantially less than the Skyrme 
value 1/3. As pointed out in ref. [19] such a reduction of the isovector spin-orbit mean- field 
(proportional to Vp„ — Vp p ) is favorable for an accurate description of isotope shifts in the Pb 
region. It seems that the isospin dependence of the spin-orbit interaction is a property for which 
two-pion exchange induced effects can play a role in nuclear structure. 

For the sake of completeness we also note that even the one-pion exchange Fock diagram 
gives rise to non-vanishing spin-orbit strength functions: 

FUkf)^ = -3G so (A: / ) (l7r) = (16 ^ )2 [^ 2 ln(l + 4 M 2 )-4] . (16) 

In order to arrive at this result one has to expand the relativistic pseudovector nNN- vertex up to 
order 1/M 2 . As expected the l7r-exchange contributions to the spin-orbit strengths are negligibly 
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small. For example, at a density of p /2 ~ 0.08 fm~ 3 one gets F^ n \kf) ~ — 0.79MeVfm 5 which 
corresponds in magnitude to less than 1% of the empirical value 90MeVfm 5 . 

Finally, we perform the same analysis for the (Vp) 2 -term in the nuclear energy density 
functional: Fy(kf) (Vp) 2 . Its strength function Fy(kf) is also decomposed into contributions 
from zero-range NN-contact interactions and long-range one- and two-pion exchange. From the 
complete set of four-nucleon contact-couplings written down in ref. [9] one gets: 



Mkf) 



(ct) 



32 



5127T 



(14Ci + C 2 - 6C 3 - 3C 4 - 2C 6 - C 7 ) 



6C{%) + 6C( 3 ^) - C( 1 P 1 ) - C( 3 P ) - 3C( 3 P 1 ) - 5C( 3 P 2 ; 



:i7) 



where we have reexpressed the relevant linear combination of C±, . . . , CV through the so-called 
spectroscopic low-energy constants. In that representation we obtain from the entries in Table 
IV of ref. [9] for the three high-precision NN-potentials CD-Bonn, Nijm-II and AV-18 the values 
F^ t] (k f ) = 82.4MeVfm 5 , 78.4MeVfm 5 and 82.5MeVfm 5 , respectively. The contribution of the 
one-pion exchange Fock-diagram to the strength function Fy(kf) reads on the other hand [8]: 

35g 2 A u- 7 f 7 21m 



Mkf) 



(lTT) 



; 18m 



(lQm 7T f 7r ) 2 
+10 arctan 2m + 
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96m ' 



16m 3 



35 

ln(l + 4m 2 ) 



}■ 



(18) 



where we have included also the relativistic l/M 2 -correction. The contributions of (lowest order) 
irreducible 27r-exchange to the NN-potential have been evaluated in ref. [11] (see eqs. (14,15) 
therein). Application of the density- matrix expansion of Negele and Vautherin [7] as outlined 
in section 2 of ref. [8] leads then to the following contribution of irreducible 27r-exchange to the 
strength function F^(kf): 



Mkf) 



(27T) 
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245(71^ + lO^ 2 - 1) + 1575m 2 (43^ + 6g 2 A - 1) 



ln 2 (M 



2u 



M z 



— ln(M + VI + v?) 
3 



%1 " I) 

xVl +m 2 [735(1 - 10^ - 7lg\) + 5m 2 (318594 + 3502^ - 1201) 
+2m 4 (869^ - 230^ - 79) + 24m 6 (2m 2 - l){llg\ - 10^ - 1) 
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Again, in order to have a clean separation between (regularization dependent) short- distance 
effects and genuine long-distance contributions we have normalized this expression to Fv(0)^ = 
0. 
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Fig. 3: The strength function Fy(kf) related to the (V p) 2 -term in the nuclear energy density 
functional versus the nucleon density p = 2k^/3n 2 . The three horizontal dashed lines show the 
constant values Fy{kf) = [9ti — (5 + 4x2)t2]/64 of the Skyrme forces Sly [17], MSk [18] and 
SIII [14]- The dashed-dotted line gives the mean value F^(k f Y c ^ = 81.1 MeVfm 5 of the three 
high-precision NN-potentials CD-Bonn, Nijm-II and AV-18. The full line includes in addition 
the In- and In-exchange contributions eqs.(18,19). 

In Fig. 3 we show the strength function F^{kf) versus the nucleon density p = 2kj/3n 2 . 
The three horizontal dashed lines represent the constant values F^(kf) = [9ti — (5 + 4a; 2 )t 2 ]/64 
of the Skyrme forces Sly [17], MSk [18] and SIII [14]. The dashed-dotted line gives the mean 
value F v (Jfe/) (ct) = 81.1MeVfm 5 of the three high precision potentials CD-Bonn, Nijm-II and 
AV-18. The full line includes in addition the In- and 27r-exchange contributions eqs. (18,19). 
One observes a rough agreement between these different curves. For densities p > 0.05 fm~ 3 our 
prediction for F^(kf) exceeds the mean value of the phenomenological Sly4-7 forces [17] only 
by about 15%. Below that density a strong rise of the strength function Fy(kf) sets however 
in. This behavior originates from the static l7r-exchange contribution with its inherent chiral 
singularity m~ 2 which becomes dominant at extremely low densities kf « m n /2 [8]. Our 
derivation of the strength function F v (/c/) is based on the density- matrix expansion of Negele 
and Vautherin [7] which has been found to become inaccurate for low and nonuniform densities 
[25]. Therefore one should not trust the full line in Fig. 3 below p < p /A ~ 0.04 fm~ 3 . It is also 
evident from Fig. 3 that the nuclear structure phenomenology does presently not constrain the 
(Vp) 2 -term as strongly than as the spin-orbit coupling term. The existing Skyrme forces show 
an appreciable variation of the associated strength parameter. 

In summary, we analyzed in this work the spin-orbit coupling term in the nuclear energy 
density functional in terms of a zero-range NN-contact interaction and finite-range contribu- 
tions from chiral two-pion exchange. We have shown that the strength of the spin-orbit contact- 
interaction as extracted from high-precision NN-potentials agrees perfectly with that of phe- 
nomenological Skyrme forces employed in non-relativistic nuclear structure calculations. The 
numerical results collected in Table IV of ref.[9] have been essential in order to establish this 
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fact. Additional long-range effects from chiral two-pion exchange seem to play a role for the 
isospin dependence of the spin-orbit coupling. The ratio of the isovector to the isoscalar spin- 
orbit strength gets significantly reduced below the Skyrme value 1/3. Such a reduction is fa- 
vorable for an accurate description of the isotope shifts in Pb nuclei [19]. Furthermore, we have 
performed the same analysis for the strength function of the (Vp) 2 -term and found values not 
far from those of phenomenological Skyrme parameterizations. 
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